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We study the phase diagram of the zero-temperature, one-dimensional Bose-Fermi-Hubbard model for fixed 
fermion density in the limit of small fermionic hopping. This model can be regarded as an instance of a disor- 
dered Bose-Hubbard model with dichotomic values of the stochastic variables. Phase boundaries between com- 
pressible, incompressible (Mott-insulating) and partially compressible phases are derived analytically within a 
generalized strong-coupling expansion and numerically using density matrix renormalization group (DMRG) 
methods. We show that first-order correlations in the partially compressible phases decay exponentially, indi- 
cating a glass-type behaviour. Fluctuations within the respective incompressible phases are determined using 
perturbation theory and are compared to DMRG results. 

PACS numbers: 



I. INTRODUCTION 

Ultracold atoms in optical lattices provide an experimen- 
tally accessible toolbox for simulating strongly correlated 
quantum systems [1-5]. The interaction of the atoms gives 
rise to local Hamiltonians on a lattice that can be character- 
ized in their microscopic details. Moreover, by means of mix- 
tures different species, Feshbach resonances, or additional op- 
tical lattices, an unprecedented control over system parame- 
ters can be achieved. Following first experiments showing a 
Mott-superfluid phase transition [3] in a bosonic system [4, 5], 
a plethora of systems of cold atoms have been studied. This 
includes mixtures of bosonic and fermionic atoms, - studied 
in experiments [9-12] and theory [13-24] - giving rise to a 
rich phase diagram and complex physics, including fermion 
pairing, phase separation, density waves, and supersolids. 

Early theoretical studies of the BFHM within mean-field 
and Gutzwiller decoupling approaches [13, 15, 16] as well as 
exact numerical diagonalization [17] revealed the existence of 
Mott-insulating (incompressible) phases with incommensu- 
rate boson filling. In comparison to the Bose-Hubbard model 
where the incompressible phases are entirely characterized by 
the local boson number, the corresponding phases for Bose- 
Fermi mixtures display a much richer internal structure. A 
rather complete description of these phases can be obtained 
using a composite-fermion picture [14], which predicts den- 
sity waves with integer filling, the formation of composite- 
fermion domains (phase separation), composite Fermi liquids 
and BCS-type pairing. The properties of ID Bose-Fermi mix- 
tures in the compressible phases were analyzed in terms of 
fermionic polarons using a bosonization approach [19]. In 
two spatial dimensions the existence of super-solid phases 
was predicted [21] which is characterized by the simultane- 
ous presence of a density wave and long-range off-diagonal 
order for the bosons. The persistence of a density wave with 
noninteger fillings in the compressible phases was shown in 
[18]. There are also a few exact numerical studies using both 
quantum Mote-Carlo [23] and DMRG calculations [24]. 

In the present work, we consider a Bose-Fermi mixture 
in a one-dimensional, deep periodic lattice described by the 
Bose-Fermi-Hubbard model (BFHM). In particular we study 
the case of small fermionic hopping, where the presence or 



absence of a fermion at a lattice site results in a dichotomic 
random alteration of the local potential for the bosons. We 
show that for this limiting case a rather accurate prediction 
of the incompressible (Mott-insulating) phases is possible us- 
ing a generalized strong-coupling approach. To verify this ap- 
proach we perform numerical simulations using the density- 
matrix renormalization group (DMRG) [27]. We predict the 
existence of partially compressible phases and provide numer- 
ical evidence that they have a Bose-glass character. Finally 
we calculate local properties in the incompressible phases and 
draw conclusions about the validity of effective theories. 



H. THE MODEL 

We consider a mixture of ultra-cold, spin polarized 
fermions and bosons in an optical lattice. In the tight binding 
limit of a deep lattice potential, the system can be described 
by the Bose-Fermi Hubbard model [13]. We here consider 
a semi-canonical model, in which the number of fermions is 
hold constant, but in which we allow for fluctuations of the 
total number of bosons determined by a chemical potential \i. 
The corresponding Hamiltonian reads 
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Here, Cj and bj are the annihilation operators of the fermions 
and bosons at lattice site j, respectively, and 
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the corresponding number operators. The par- 
ticles can tunnel from one lattice site to a neighboring one, 
the rate of which is described by Jg and Jp for bosons and 
fermions, respectively. V is the on-site interaction strength 
between the two species, while U accounts for intra-species 
repulsion of bosons, which will define our energy scale and 
we set henceforth U = 1. 
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Throughout the present work, we will focus on the case of 
heavy, immobile fermions, i.e., we consider the limit in which 
Jf = is a good approximation. In this case the effect of 
the fermions reduces to a dichotomic random potential at site 
j for the bosons, depending on whether a fermion is at site j 
or not. This means that the local potential is altered by 



V, if a fermion is present at site j, 
0, otherwise. 



(2) 



We will systematically investigate to what extent this limit of 
the Bose-Fermi-Hubbard model can be described as an spe- 
cific instance of a disordered Bose-Hubbard model, 
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We will see that this simple model shows on the one hand im- 
portant features of the full Bose-Fermi-Hubbard model. On 
the other hand, we will see that this leads to important quali- 
tative differences to the phase diagram of the disordered Bose- 
Hubbard model with continuously distributed on-site disorder, 
as studied in Refs. [4, 8]. Depending on the physical situa- 
tion of interest we will consider two cases of disorder: If the 
fermionic tunneling is small but sufficiently large such that 
on the time scales of interest relaxation to the state of total 
minimum energy is possible, the fermion induced disorder is 
referred to as being annealed. In this case the ground state is 
determined by minimization over all possible fermion distri- 
butions. If the fermion tunneling is too slow or the temper- 
ature too high the disorder is an actually random distribution 
called quenched. 



III. COMPRESSIBLE AND INCOMPRESSIBLE PHASES 

In this section we derive the phase diagram of the BFHM 
with immobile fermions. More specifically, we will approx- 
imate the boundaries between compressible and incompress- 
ible phases employing a generalization of the familiar strong- 
coupling expansion [25] to the present case of bosons with 
a modified potential due to the presence of fermions. This 
will be compared to the predictions of several instances of 
mean-field approaches [15, 16]. Furthermore, a comparison 
with numerical results in one spatial dimension obtained by a 
DMRG computation will be given. Our strong-coupling ex- 
pansion reveals the existence of novel phases whose character 
will be discussed in the subsequent section. 



A. Ultra-deep lattices 

We first discuss the simple case of an ultra-deep lat- 
tice for the bosons, such that their hopping can be ne- 
glected. In this situation where Jp = Jb = is a 
good approximation, the Hamiltonian becomes diagonal in 



the occupation number basis. This basis will be denoted as 
{|ni, • • • , rijv")|mi, • • • , mjv)}, where m,j = 0, 1 denotes the 
number of fermions at site j and rij = 0, 1, ... the correspond- 
ing number of bosons. N labels the total number of lattice 
sites in the one-dimensional system. The problem of find- 
ing the ground state reduces to identifying product states with 
the lowest energy. By fixing the total number of fermions 
Np = Ngp, this amounts to minimizing 

E = \ Y n i( n j - !) - (m - v) y n 3 - V Y n i- 



where T denotes the set of qfN = Np sites with a fermion 
and AT the set of (1 — qf)N = N — Nf sites without a 
fermion, g p denoting the fermionic filling factor. The energy 
is obviously degenerate for all fermion distributions and the 
ground state is given by an equal mixture of all states with 
state vectors 



|V>0> 



K,i> 



K,o). 



(4) 



Here, 



m = max 



{0, [1/2 + 0* - V)]j, n = max{0, [1/2 + fj]j, 

(5) 

is the local boson number for sites with (J 7 ) or without (AO a 
fermion and [.] denotes the closest integer bracket. In other 
words, the degenerate states with lowest energy will have 
QfN sites with ni bosons and one fermion and N(l — qf) 
sites with ?i bosons and no fermion. For the case of zero or 
unity fermion filling, Qp = 1 the situation becomes particu- 
larly simple as we encounter the pure Bose-Hubbard model 
with an effective chemical potential u cti = fx — Vqf- 

Since m and no are integers there are adjacent intervals 
of [i where the occupation numbers do not change. In these 
intervals the system is incompressible, i.e., 



d(Ej nj) 
dfj, 



= 



(6) 



and the points between two intervals are quantum critical 
points. This behavior, illustrated in Fig. 1, is very similar to 
that of the Bose-Hubbard model except that here the bosons 
can be incompressible even for non-integer filling qb as we 
have qb = no + £>f(«i — no). Following Ref. [14] we la- 
bel the difference no — n\ in the bosonic number mediated 
through the presence of a fermion by s. The local ground state 
can either consist of no bosons and no fermion or n\ = no — s 
bosons and one fermion. These state vectors will be denoted 
as|no,0) = |0) and \tiq — s, 1) = The value of s de- 
pends on fx and V and can be a positive or negative integer. 
Both these vectors are eigenvectors of the number operator 



h = n i + sm i 



(7) 



with the same integer eigenvalue no and (AQ|) = 0. Thus 

incompressible phases have a commensurate number Q and 
can be characterized by the two integers no and s Since no 
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FIG. 1: (Color online) Phases of BFHM for Jb = Jf = for differ- 
ent inter-species couplings < V < 1 (lowest diagram), 1 < V < 2 
(middle diagram), — 1 < V < (upper diagram), and U = 1. n in- 
dicates the number of bosons (empty circles) at the site, m the num- 
ber of fermions (red filled circles). The horizontal red bars illustrate 
the boson number ni for sites with a fermion (m = 1) as function 
of the chemical potential, the horizontal blue bars correspondingly 
the boson number no for sites without a fermion (m = 0), which 
is identical to the BHM. The values of fi where a transition between 
different boson numbers no occurs either at sites without a fermion 
(m = 0) or with a fermion (m = 1) are quantum critical points. 



and ri\ are integers and increase monotonically with [i , there 
is a jump in the total number of bosons when moving from 
one incompressible to the adjacent one. All systems with bo- 
son number in between these values are critical and have the 
same chemical potential since Jb = 0. The average boson 
number per site in the incompressible phases does not have 
to be an integer, however. The existence of Mott phases with 
non-commensurate boson number is a direct consequence of 
the dichotomic character of the fermion induced disorder. A 
similar behavior has been predicted for superlattices, which 
can be considered as dichotomic disorder in the special case 
of anti-clustering [28, 29]. In general Mott-insulating phases 
with incommensurate boson numbers exist for any disorder 
distribution that is non-continuous. 



B. Minimum energy distribution of fermions for small bosonic 
hopping 

In order to understand the physics for disorder due to the 
presence of fermions we need to discuss the influence of the 
distribution of fermions to the ground state energy. The ener- 
getic degeneracy of different fermion distributions in the in- 
compressible phases is lifted if a small bosonic hopping Jb is 
taken into account. Near the quantum critical points the boson 
hopping leads to the formation of possibly critical phases with 
growing extent. We first restrict ourselves to regions where in- 
compressibility is maintained, i.e. sufficiently far away from 
the critical points. 
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FIG. 2: (Color online) Fermion distributions ordered by ground state 
energy. Blue - lowest energy, red -highest energy for Jf — 0, Jb = 

0. 02, U = 1. top: attractive boundary, (i) V = 1.5, no = 1, s = 1 

1. e. K cS = -0.002, and (ii) V = 1.5, n = 2, s = 1 i.e. K cS = 

0. 001. bottom: repulsive boundary (iii) V = —1.5, no = 0, s = —1 

1. e. K eS = -0.002, (iv) V = -1.5, n = l,s = -1 i.e. K cfi = 
0.001. 



In order to obtain a qualitative understanding of the effects 
of a finite bosonic hopping we have performed a numerical 
perturbation calculation on a small lattice. Fig. 2 shows differ- 
ent distributions of 4 fermions over a lattice of 8 sites ordered 
according to their energy for different parameters in 6th order 
perturbation. 

One notices that the lowest energy states are either given by 
fermion distributions with maximum mutual distance (anti- 
clustered configuration) or minimum mutual distance (clus- 
tered configuration) modified by boundary effects. This be- 
havior can in part be explained by the composite fermion pic- 
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ture introduced in [14]. The composite fermions are defined 
for the phase (n , s) by the annihilation operators: 



fi 



fi = 



(n -s)\ 
(no - s)\ 



C i j 



for s > 0, 



for s < 0. 



(8) 



(9) 



For each no and s = 1, the full BFH-Hamiltonian, Eq. (1), 
with Jp = gives in second order in Jb rise to the effective 
Hamiltonian [14] 

U-k) 

where (., .) denotes nearest neighbors. Here, as Jp = 0, we 
find the effective coupling (note that again, U = 1) 

■n (n + l-s) (n -s)(n + l) 



1-s + V 1 + s-V 

-n (n + 1) - (n - s)(n + 1 - a) 



(11) 



Composite fermions cannot occupy the same lattice site, but 
there may be nearest neighbor attraction (K e g < 0) or repul- 
sion (K e g > 0). Associating a site with a composite fermion 
with a spin-up state and a site without a fermion with spin 
down, Eq. (10) corresponds to the classical Ising model with 
fixed magnetization and anti-ferromagnetic (K c s > 0) or fer- 
romagnetic coupling {K c ff < 0). 

As a consequence, to this order in perturbation theory, if 
Kefi < 0, the energy is smallest for fermion distributions that 
minimize the surface area of sites with and without a fermion 
(referred to as clustering). In this setting, we can take the 
fermion distribution to form a block of occupied sites. 

The other regime is the one for K c g > 0. Then, the 
fermions repel each other, and they form a pattern with maxi- 
mum number of boundaries for small Jb, referred to as anti- 
clustering. That the fermions attain a distribution with max- 
imum distance cannot be explained by the effective model 
due to its perturbative nature. In all of our numerical simula- 
tions using the density matrix renormalization group (DMRG) 
we found however that a positive K c g always lead to anti- 
clustering with maximum distance. 

The ground state energies of the various fermionic distribu- 
tion differ only by a small amount which is on the order of 
J%/U or even higher powers. Also for temperatures which 
are still small enough to treat the bosonic system with given 
disorder as an effective T = problem, but larger than the 
energy gap between different fermion distributions, i.e. for 
Jb(Jb /U) n <C ksT <C Jb, the various fermion distribu- 
tions will be equally populated. Thus it seems more natural 
to consider the case of quenched, random disorder rather than 
that of annealed disorder. 



C. Compressible and incompressible phases for finite Jb 

We now discuss the boundaries of the incompressible 
phases for finite bosonic hopping. To this end we extend the 



strong coupling expansion of Ref. [25] and complement the 
results with numerical DMRG simulations. The strong cou- 
pling expansion provides a rather accurate description for the 
Bose-Hubbard model even on a quantitative level. 

Let us consider a phase with (no,s) and Np = QpN 
fermions, i.e., a phase with Np sites containing no — s bosons 
and a fermion and N — Np sites with no bosons. The ground 
state vector for Jb = is then found to be 




- f\ (no— s) 

|V'o) 

The energy density is given by 

U 

£0 = 2 



,0)|0,.. 



,0). 
(12) 



(1 - QF)no{n 
Vg F (n - s). 



1) + QF(n - s)(no - s - 1) 

(13) 



We now consider states with a single additional boson 
(bosonic hole). In contrast to the actual Bose-Hubbard model 
in the absence of fermions, we here have to distinguish two 
cases, where a boson (bosonic hole) is added to a site with a 
fermion. Up to normalization, we have 

\^+^y = Oj\tpo), i^-^y = o,j\ipo), j e (i4) 

or without a fermion 

\^+mY = a]\ipo), \4>-.Af) J = o,j\ipo), jeJV. (15) 

All of these vectors are eigenvectors of the BFH-Hamiltonian 
for Jb = with respective energies 

= E + V + U(n -s), (16) 

E-s = E -V + U(n -s-l), (17) 

E +M = E + Un , (18) 

E_.x - E + U(n -1), (19) 



where Eq = Neq. 
read 



The corresponding chemical potentials 



E+.jr -E = V + U{n - s), (20) 
= E - E_ t r = - U, (21) 



and 



— E + _f S f - E 
M-jv = E — E-^ 



Un , 



(22) 
(23) 



Except from the special case V = Us, the energies E± f 
and E± t tf all differ from each other. Thus we can determine 
the phase boundaries for Jb ^ by degenerate perturbation 
theory within the subspaces j E T and j E AT separately. 

There will be a second order contribution in Jb for sites 
j that have at least one neighboring site of the same type. 
For isolated sites degenerate perturbation theory will lead 
only to higher order terms in 0(J B ). Since the boundaries 
of the incompressible phases are determined by the overall 
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FIG. 3: Color online) Phase diagram from strong-coupling expansion 
and U = 1,V — 1.5. Red areas (A) indicate truly incompressible 
Mott regions with gaped particle-hole excitations everywhere. Green 

(B) or blue (C) areas are partially compressible quasi-Mott regions 
with gaped particle-hole excitation for sites with (without) a fermion 
but ungapped excitation in the complementary region. 

lowest-energy particle-hole excitations, we can construct the 
expected phase diagram in the case of extended connected re- 
gions of fermion sites coexisting with extended connected re- 
gions of non-fermion sites. In this case we can directly apply 
the results of Ref. [25] to sites with and without fermions 

M±,:f/AT = /4 ) jF/Ar + (5 M±( ri o, J b) (24) 

where 

Sp + (n ,J B ) = -2J B (n + 1) + J B n\ 

+ j|n (n + l)(n + 2), (25) 

5fi-(n , Jb) = 2J B n Q - J|(n + l) 2 

- 4no(ng-l). (26) 

This gives rise to two overlapping sequences of quasi-Mott 
lobes shifted by the boson-fermion interaction V as shown in 
Fig. 3. 

The system is truly incompressible only in the overlap re- 
gion of the quasi-Mott lobes (A). Points which are within one 
of the two sequences of quasi-Mott lobes but not in both (cases 
B or C) are partially incompressible with an energy gap for a 
bosonic particle-hole excitation on a site with (B) (without 

(C) ) a fermion but without a gap for a corresponding excita- 
tion on a complementary site. The properties of these partially 
incompressible phase will be discussed later. 

These strong coupling results will now be complemented 
by numerical calculations using a DMRG simulation for a 
system with fixed fermion distribution an open boundary con- 
ditions. The local Hilbert space for the bosonic sector is 
span{|0), . . . , |6)}, so it is truncated at 6 bosons. The DMRG 
computation is done for both clustered and anti-clustered 
fermion distributions. The corresponding graphs for the phase 
boundaries are shown in Fig. 4. 

One recognizes nearly perfect agreement between numerics 
and strong-coupling prediction in the case of clustering. This 
is expected since in the clustered case the majority of sites 




bosonic hopping J 



FIG. 4: (Color online) Comparison of strong-coupling approxima- 
tion (full line -) and DMRG for boundaries of incompressible phases 
for fixed distribution of fermions corresponding to clustering (in- 
ner crosses, X) or anti-clustering with maximum distance (outer 
crosses ,x) . V = 1.5 (top figure) and V = —1.5 (bottom figure). 
qf = 0.25, and U = 1. 



has neighbors of the same type. In the case of anti-clustering, 
however, the incompressible lobes extend much further into 
the region of large boson hopping with a critical J B of about 
1 for a fermion filling of qf — 1/4 at 1^ = —1.5. The latter is 
to be expected since in this case hopping to nearest neighbors 
is suppressed if the neighboring sites are of a different type 
(J- or NT). Here the curves of the critical chemical potential 
Merit ( Jb ) that correspond to a bosonic particle-hole excitation 
at a fermion site (here /i C rit(0) = —1.5,-0.5,0.5, 1.5 etc.) 
start with a power Jg determined by the minimum number of 
hops required to reach the next fermion site, i.e. 7 = 1/pf, 
if Qf < 1/2. If the fermion filling is larger than 1/2 the 
picture changes and the non-fermion sites (hole sites) cause 
Mcrit(^s) ~ J B with 7 = 1/(1 — Qp). In principle it is 
possible to extend the strong-coupling perturbation expansion 
to any fermion distribution, which is however involved. Fig. 
5 shows the prediction of a cell-strong coupling expansion 
[30] for an anti-clustered, fixed fermion distribution which is 
equivalent to bosons in a super-lattice potential [36]. 
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FIG. 5: (Color online) Comparison of cell-strong coupling approxi- 
mation (full line -) obtained from [30] and DMRG for boundaries 
of incompressible phases for fixed distribution of fermions corre- 
sponding to anti-clustering with maximum distance (full dots,*) . 

V = -1.5, qf = 0.25, and U = 1. 

We now want to argue that the strong-coupling expansion 
for a clustered fermion distribution provides an accurate pre- 
diction for the boundaries of the incompressible phases in the 
case of quenched, random fermion disorder. Since in the ther- 
modynamic limit any local distribution of fermions is real- 
ized at some places in the lattice, the actual phase boundaries 
are determined by the fermion configuration that leads to the 
smallest incompressible regions. Since this is the case for 
a clustered fermion configuration, which in turn is well de- 
scribed by the strong-coupling expansion, the latter gives a 
rather accurate description of the phase transition points be- 
tween compressible and incompressible phases. 

For the case of annealed fermion distribution the strong- 
coupling expansion is expected to give only less accurate re- 
sults. This can be seen from Fig. 6 where we compare the 
predictions of the strong-coupling approximation with those 
from a DMRG simulation for annealed fermionic disorder and 
a mean-field ansatz. Within the mean-field approach, e.g. of 
[15], hopping is included to the system as a perturbation to the 
ground state 

Is > = V 1 - 0fK,o) + VeFK (27) 

Using this ground state and introducing a global bosonic or- 
der parameter ip, the phase boundaries can be found using 
the usual Landau argumentation. For details see [15, 16]. 
Fig. 6 shows the resulting phase diagram compared to DMRG 
data for annealed disorder and strong-coupling predictions. 
When comparing the different data sets one recognizes that 
the mean-field predictions are qualitatively correct but as ex- 
pected only moderately precise quantitatively. It should be 
mentioned that the accuracy of the mean-field approach be- 
comes worse even for Jb — > for a disorder with maximum 
anti-clustering. The numerical data were obtained by letting 
the DMRG code freely evolve in the manifold of fermionic 
distributions. The obtained distribution then gives a state 
which is at least close to the ground state. Since this pro- 
cedure is prone to get stuck in local minima we checked the 
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FIG. 6: (Color online) Incompressible phases for an annealed 
fermion distribution and U = 1, V = 0.25, qf = 0.25. Shown 
are the strong-coupling results (outer full line -), the mean-field re- 
sults from [15] (inner full line -) and results from a DMRG calcula- 
tion with Jf = (crosses x). Vertical line indicates position of the 
density cut shown in Fig. 9 in the following section. 



consistency of our results by implementing different sweep 
algorithms. In these algorithms the fermionic hopping was 
not taken to be zero but was given a finite initial value which 
was decreased during the DMRG sweeps to the final value 
zero. To ensure proper convergence we compared the data for 
a few representative points (Jb = 0.07 boundaries of (no = 
l,s = 1) lobe; Jb = 0.15 boundaries of (no = l,s = 0) 
lobe; Jb = 0.03 boundaries of (no = 2, s = 0) lobe) to the 
data obtained from two different sweep strategies [37]. The 
difference in the chemical potential is of the order of 3% in- 
dependent of the sweep strategy an therefore negligible on the 
scale of the plot. 



D. Influence of finite fermionic hopping 

The question arises how the phase diagram changes if a 
finite but small fermionic hopping is included. The case 
Jf 7^ should be compared to the case Jp = for annealed 
fermionic disorder. Fig. 7 shows a comparison of DMRG data 
for Jp = and Jp = Jb- One recognizes that the influence 
of a small fermionic hopping is rather small. 



E. Finite size extrapolation 

The DMRG simulations are done for finite lattices and thus 
finite-size effects influence the results. To eliminate these ef- 
fects each data point is obtained by a finite size extrapola- 
tion. This is particularly important if one wants to determine 
the critical values of Jb for the compressible-incompressible 
transition. Figure 8 shows the extrapolation of the tip of the 
lowest Mott phase in Fig. 7 for Jp = Jb to infinite lattice 
sizes N — > oo. From a fit of J c to log(iV) we find the criti- 
cal point in the thermodynamic limit J c = 0.16038. The data 
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FIG. 7: (Color online) Incompressible regions for U = 1, V = 0.25, 
qf = 0.25. Shown are the results from a DMRG calculation with 
Jf = and annealed fermionic disorder (X ) and Jf = Jb (° )• The 
phases are the same as in fig. 6. Vertical line indicates position of the 
density cut shown in Fig. 12. 
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FIG. 8: (Color online) Thermodynamic limit extrapolation for the 
critical point of the no — 1, s = 1 lobe (qf = 1/4, qb = 3/4) in 
figure 7. The critical point is found at J c = 0.16038 



for different system lengths show the expected 1 /N behavior 
shown in Ref. [26] for the BHM. 



IV. PARTIALLY INCOMPRESSIBLE PHASES 

A. Limit of vanishing fermionic hopping 

Within the strong-coupling approximation discussed in the 
previous section we have identified regions in the (j,—Jb phase 
diagram where bosonic particle-hole excitations are gapless if 
they occur on a fermion (non-fermion) site but have a finite 
gap on a complementary i.e. a non-fermion (fermion) site. 
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FIG. 9: (Color online) Density cut along the vertical line in figure 6. 
From top to bottom: overall average density (-), average density for 
sites without a fermion (-), average density for sites with a fermion 
(-). Inset: Dependence of particle number at one particular site as 
function of chemical potential without averaging. 



Associated with this is a partial incompressibility 



0. 



0, 



or 



7^0, 



(28) 



^0. 



This is illustrated in Fig. 9. Here the average boson number 
per site obtained from a DMRG simulation with annealed dis- 
order is shown as a function of the chemical potential for 
constant bosonic hopping. The curve corresponds to the pa- 
rameters of Fig. 6 for the vertical cut shown in that figure at 
Jb = 0.02. Also shown are the corresponding values only 
for fermion sites and non-fermion sites respectively. In the 
partially compressible phases the average boson number in- 
creases only for one type of sites while it stays constant for 
the other. In the DMRG code the energy per particle is calcu- 
lated as a function of the total number of bosons N which then 
yields the chemical potentials fx + (N) = E(N + 1) - E(N) 
and H-{N) = E(N) — E(N — 1). Averaging over few val- 
ues of N in the compressible phase is needed here since the 
ground state fermion distribution changes with changing bo- 
son number leading to a non-monotonous dependence of fj, on 
the boson number. 

We now discuss the properties of the single-particle den- 
sity matrix (oja; +m ) in the partially incompressible phases. 
For very large values of Jb the system is expected to have 
a Luttinger-liquid behavior in Id and to possess long-range 
off-diagonal order in higher dimensions. In Id we expect 
that the Luttinger-liquid behavior disappears in the partially 
incompressible phases and that correlations decay exponen- 
tially. This is because in this case a single (static) impurity is 
sufficient to prevent the build-up of long-range correlations. 
In higher dimensions there will be a critical fermion (or hole) 
filling fraction above which off-diagonal order is suppressed. 
This critical fraction is determined by percolation thresholds 
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FIG. 10: (Color online) DMRG simulations (X) of first-order cor- 
relations (a\ai +m ) for V = 0.25, qf = 0.25 and (7 = 1 for a 
lattice of 512 sites and Nb = 448 bosons in the case of annealed 
disorder, (top curve:) Jg = 0.07, line corresponds to exponential fit 
oc exp{— m/lc} with l c = 1.7. The exponential decay for small Jb 
is apparent, (bottom curve:) Jb = 0.2, line corresponds to algebraic 
fit oc m~ a with exponent a = 0.33. 



and for annealed fermionic disorder depends on the actual 
fermion distribution in the ground state (e.g. clustered or anti- 
clustered). For a random fermion distribution in 2d the thresh- 
old is gf [t = 0.5927 (or 1 - gf n = 0.5927 if non-fermion 
sites are incompressible). The corresponding number for 3d 
is g c f lt = 0.3116. 

Fig. 10 shows the first-order correlations (ajai+ TO ) as func- 
tion of the distance m for an annealed fermion distribution 
obtained from DMRG simulations for a rather large lattice of 
512 sites with incommensurate boson filling (Nb = 448) and 
qf = 1/4. For Jb = 0.07 strong exponential decay with cor- 
relation length l c = 1 .7 is found corresponding to a glass-type 
behavior, while for Jb = 0.2 correlations decay algebraically 
with m -0 ' 33 , which corresponds to a Luttinger liquid. Note 
that for the chosen boson number, which corresponds to a non- 
commensurate value of Q there is no incompressible phase. 

Fig. 11 shows the first-order correlations for a random, 
quenched fermion distribution averaged over 100 realizations 
with non-commensurate boson number (qb = Nb/N = 



FIG. 11: (Color online) DMRG simulations (x) of first-order corre- 
lations (a\a,i+m) for V = 1.5, qf — 0.375 and (7 = 1 for a lattice 
of 128 sites and Nb = 184 bosons averaged over 100 fermion distri- 
butions, (top curve:) Jb = 0.03, line corresponds to exponential fit 
oc exp{— m/lc} with l c = 2.9. The exponential decay for small Jb 
is apparent, (bottom curve:) Jb = 0.2, line corresponds to algebraic 
fit oc m~ a with exponent a = 0.37. To avoid finite size effects at 
short and long ranges, only the sites in between 13 and 110 are taken 
into account. 



184/128). Despite the sampling noise one recognizes the 
transition between exponential decay with correlation length 
l c = 2.9 for Jb = 0.03, and a power-law decay with to -0,37 
for Jb = 0.2 corresponding to a Luttinger liquid. Jb = 0.03 
is within a partially incompressible phase, Jb = 0.2 outside. 

The numerical results and the above discussion indicate that 
the partially incompressible phases have a glass-type charac- 
ter. A detailed discussion of the Bose-glass to superfluid tran- 
sition will be given elsewhere [33]. 



B. Small fermionic hopping 

If there is a non-vanishing but small fermionic hopping, 
partial incompressibility is lost. Still the increase of the bo- 
son number with increasing chemical potential at one type of 
sites is substantially less that that on the complementary type 
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FIG. 13: Possible single-hop excitations; the energies are given by 
A x i? to A 4 E (see text for definitions). Filled circles are fermions, 
open circles bosons. 



FIG. 12: (Color online) Density cut for the same parameter as Fig.9 
but for Jb — Jf = 0.05 (see Fig.7). 



of sites. 



(9// 



9// 



or 



(29) 



Fig. 12 shows the density cut obtained from DMRG simula- 
tions for the parameters of Fig.9 but for Jb = Jf = 0.05. It 
should be noted that in contrast to Fig.9 averaging over sites 
is not needed due to the finite mobility of the fermions. The 
simulations show that the glass-type character of the phases 
survives. We expect a crossover from glass-type to Luttinger 
liquid behavior with increasing fermionic hopping. In addi- 
tion due to the stronger back-action of the boson distribution 
to the fermion distribution other phases such as density waves 
emerge [18]. A discussion of the Bose-Fermi Hubbard model 
in the limit of large fermion mobility will be given elsewhere 
[35]. 



V. FLUCTUATIONS 

In this section, we will determine the fluctuations of the 
bosonic number operator for vanishing fermionic hopping 
Jf = inside the quasi Mott lobes for quenched disorder. 
To this end, second order perturbation theory will be applied 
and compared to numerical results from DMRG. 

For vanishing bosonic hopping the ground state with a fixed 
number of fermions is clearly highly degenerate, from the dis- 
tribution of Np fermions in a lattice with N sites. For the case 
of quenched disorder with fixed positions of fermions, consid- 
ered here, this degeneracy is inconsequential. This allows to 
develop a tractable approach based on non-degenerate pertur- 
bation theory for a given fermion distribution and subsequent 
averaging over all of these distributions. In order to evalu- 
ate the fluctuations of the bosonic number operator, we hence 
have to determine 



which is independent of the lattice site j due to translational 
invariance. Here, the classical average E is taken with respect 
to the fermionic distributions, so the average over the (^ ) 
different distributions with equal weight. 

We can hence proceed as in Ref. [34] to compute the fluc- 
tuations in the boson number, for each fermion distribution, 
followed by the appropriate average. In second order pertur- 
bation theory in Jb at Jf = 0, only bosonic hoppings to 
nearest neighbors contribute. On such two sites, clearly, four 
different situations can arise, dependent on whether or not a 
fermion is present at each of the two sites, see Fig. 13. The 
change in energy due to these excitations is given by (here, we 
have no longer taken U — 1) 



A^E = A^E = -U, 
A {2) E = -U(l -s) - V, 
A (3) £ = —U(l + s) + V, 



(31) 
(32) 
(33) 



where the superscript denotes the type of process according to 
Fig. 13. With this we are now able to calculate the fluctuations 
E(An|) of the bosonic number operator. After a number of 
steps, following the procedure of Ref. [34], we find 

E(Anf) = n {n Q + 1)(1 - g F ) 2 

Jb 
U 



+ 2z — (n - s)(n - s + l)gp 



2zJ 



2 n (n - s + 1) + (n - s)(n + 1) 

B 



U 2 - (Us - Vf 



x Qf(1~Qf), 



(34) 



(30) 



where z gives the number of nearest neighbors. The fluctua- 
tions show the expected quadratic dependence on the hopping 
strength. Moreover, in the two limiting cases gp = and 
Qp = 1 this expression coincides with the pure BHM result 
from Ref. [34]. Fig. 14 shows the analytical result compared 
with DMRG calculations for annealed disorder. For small Jb 
the agreement is rather good with increasing disagreement for 
bigger Jb, where second order perturbation theory starts to 
fail. 

Fig. 15 shows the dependence of the fluctuations of the 
fermionic density pp at a fixed hopping Jb- Also shown is 
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bosonic hopping J g 

FIG. 14: (Color online) Bosonic number fluctuations in the Mott- 
insulating lobes for different lobes at fixed density for annealed dis- 
order. The upper three lines are for V = 1.25, qf = 0.375, s — 1 
with no = 1(>), 20) and 3(>) scaled by 10 3 . Lower three lines: 
V = 1.7, g F = 7/16, s = 2 with n = 2(«), 3(») and 4(»); solid 
lines are the corresponding analytic curves 




0.3 0.4 0.5 0.6 

fermionic density p B 

FIG. 15: (Color online) Bosonic number fluctuations in the quasi 
Mott lobes for different lobes at fixed hopping. The upper three lines 
are the results for V = 1.25, J B = 0.004, s = 1 with n = 1 (►), 2 
(►) and 3 (►), shifted by 10~ 3 upwards; lower two lines: V = 1.7, 
Jb = 0.004, s — 2 with no = 2 (•) and 3 (•); solid lines are the 
corresponding analytic curves. Red circles (o) are the corresponding 
fluctuations for a clustered disorder (only for the uppermost plot). 



one numerical curve obtained with annealed disorder. The 
agreement between the analytical expression (34) and the nu- 
merical data shows, that the above derivation gives a good 
estimate for the fluctuations in the system for small bosonic 
hopping. 



the one-site energy. An important difference to the disordered 
Bose-Hubbard model [4] lies however in the distribution of 
on-site energies which is here not continuous but binary, cor- 
responding to the presence or absence of a fermion at a given 
site. As a consequence there are no extended compressible 
phases for vanishing bosonic hopping. Instead incompress- 
ible phases with in general incommensurate boson number 
emerge similar to the case of a super-lattice [28, 29]. These 
Mott-insulating phases which can be characterized by two in- 
teger parameter no and s, denoting the number of bosons at 
sites without a fermion and the shift of this number due to 
the presence of a fermion have been predicted before within 
mean-field and Gutzwiller approaches [14, 16]. Here we de- 
termined the extend of these phases using a modified strong- 
coupling expansion and numerical simulations employing the 
density-matrix renormalization group (DMRG). We showed 
that the shape of the quasi-Mott lobes depends on the actual 
fermion distribution. The latter is determined by the prepa- 
ration technique. If the fermionic hopping is small but suffi- 
ciently large such that the fermions have time to find the ener- 
getically lowest configuration, one has an annealed fermionic 
disorder, otherwise the distribution is random and frozen. For 
the annealed case we showed that in the limit of small, but 
nonzero bosonic hopping, Jb, the fermions form either a clus- 
tered or an anti-clustered configuration with maximum mu- 
tual distance. A partial explanation for this behavior could 
be found in terms of the composite-fermion model of [14]. 
For the case of random, quenched fermion distributions we 
could derive semi-analytic predictions for the phase bound- 
aries of the incompressible phase using a strong-coupling ap- 
proach [25] which agreed very well with numerical simula- 
tions. Within this approach we also identified partially com- 
pressible phases where particle-hole excitations at one type of 
site, i.e. either with or without a fermion, are gap-less, while 
the corresponding excitations at the complementary type of 
sites are gapped. The partial compressibility of these phases 
was verified by numerical simulations. We also showed that 
that the presence of partial compressibility lead to Bose-glass 
phases, which are gap-less but for which first-order correla- 
tions decay exponentially. We discussed the influence of a 
finite bosonic hopping on local properties in the quasi-Mott 
phases using a perturbative approach supplemented by numer- 
ical DMRG simulations. Finally we also discussed the influ- 
ence of a finite fermionic hopping. The numerical simulations 
indicate that many predictions remain valid for finite values 
of Jp even as large as Jb- A more detailed discussion of the 
limit of large fermionic hopping and the associated new phe- 
nomena such as density waves etc. will be given elsewhere 
[35]. 



VI. SUMMARY 
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